Lower branch coherent states in shear flows: transition and control 
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Lower branch coherent states in plane Couette flow have an asymptotic structure that consists of 
O(l) streaks, 0{R~ 1 ) streamwise rolls and a weak sinusoidal wave that develops a critical layer, for 
large Reynolds number R. Higher harmonics become negligible. These unstable lower branch states 
appear to have a single unstable eigenvalue at all Reynolds numbers. These results suggest that 
the lower branch coherent states control transition to turbulence and that they may be promising 
targets for new turbulence prevention strategies. 



Recent experiments indicate that the smallest ampli- 
tude necessary to trigger transition to turbulence in pipe 
flow scales with the inverse of the Reynolds number R, at 
least for a class of large scale perturbations [H, Q ■ That 
R^ 1 scaling, and other characteristics of the perturba- 
tions, are shown here to be consistent with a class of un- 
stable 3D traveling wave solutions of the Navier-Stokes 
equation recently discovered in all canonical shear flows 
[3 0, H, 0, Sj. These new coherent solutions arise 
through saddle- node bifurcations at R = R sn . At that 
onset Reynolds number, the solutions capture the form 
and length scales of the coherent structures that have 
long been observed in the near wall region of turbulent 
shear flows For R > R sn , the solutions separate into 
upper and lower branches. For relatively low R > R sn , 
a single traveling wave upper branch may capture the 
key statistics of turbulent shear flows remarkably well 
@, Here it is shown that the lower branch so- 

lutions in plane Couette flow obey the i?" 1 scaling and 
evidence is provided that these states form the 'backbone' 
of the phase space boundary separating the basin of at- 
traction of the laminar flow from that of the turbulent 
flow, and are therefore directly connected with transition 
to turbulence EK]. 

Incompressible fluid flow is governed by the Navier- 
Stokes equations 



d t v 



Vv + Vp = R~ 1 V 2 v, 



0, (1) 



where v(r,t) is the fluid velocity at point r and time 
f > 0, p(r, t) is the mechanical pressure that enforces in- 
compressibility and R > is the Reynolds number which 
is a non-dimcnsionalizcd inverse viscosity. The mean flow 
is in the e x direction in a channel with parallel walls at 
y = ±1. Plane Couette flow (PCF) is driven by the mo- 
tion of these walls so v = ±e x at y = ±1, for all x, z, t, in 
which case v = ye x is the laminar solution of ([T]). That 
solution is linearly stable for all R > 12]. Periodic 
boundary conditions are imposed in the wall-parallel di- 
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FIG. 1: Amplitude of x-Fourier modes for a 3D steady state 
in plane Couette flow vs. R for (a, 7) = (1,2). Top to bot- 
tom: O(l) streak uo(y, z)— u(y), O(R~ ' 9 ) fundamental mode 
|u>i|, 0(R~ 1 ) streamwise rolls (vo,wo) and o(R~ 1 ) \v2\ and 
I V3 1 . Continued beyond R — 6168 by dropping all harmonics. 
Rsn ~ 164 is the turning point where lower and upper branch 
solutions coalesce. 



rections x and z with fundamental wavenumbers a and 
7, respectively. Technical details can be found in Q. 

For traveling wave solutions, the velocity field is 
Fourier decomposed in the ^-direction as 



v(r,t) = v (y,z) 




v n (y,z) + c.c. 



(2) 



where 9 = a{x — ct), c is the constant wave velocity and 
c.c. denotes complex conjugate. The 0-mode Vo(y,z) = 
(uq, vo, Wo) consists of streamwise rolls (0, vq,wq) with 
d y v + d z w = kinematically decoupled from the 
streamwise component uq. The latter consists of an x and 
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FIG. 2: Contours of vo(y, z) (solid, top) and z)\ (solid, 

bottom) both with contours of uo(y,z) = [—2 : 2]/3 (dashed) 
for (q,7,_R) = (1,2,50171). The critical layer uo(y,z) = is 
shown as a bold solid curve in both plots. 



z averaged mean flow u(y) and streaks uo(y, z) — u(y). 

Symmetric lower branch traveling waves in plane Cou- 
ette flow (for which c = 0) have been continued to high 
R by Newton's method as in 0,0, E|- Figure Q] shows the 
scaling of the amplitudes of the various elements con- 
stituting such solutions as functions of R. The streaks 
uo(y, z) —u(y) tend to a non-zero constant while the am- 
plitude of the rolls (0, vo(y, z),wo(y, z)) scales like 1/R as 
R — * oo. The fundamental mode v%(y, z) has an approx- 
imate R~ 9 scaling, while the 2nd and 3rd harmonics 
scale approximately like R~ 16 and R~ 22 respectively. 
Higher harmonics decay faster and are not shown. This 
separation between the harmonics suggests that the 2nd 
and higher harmonics become insignificant for large R. 
Indeed, the solution was continued beyond R = 6168 
by dropping all harmonics with no significant change 
(none detectable on fig. [I]). This is unusual: as R is in- 
creased, the numerical resolution can be decreased in the 
x-direction. This is only true for the continuation of the 
lower branch solutions, and the catch is that the struc- 
ture of the lower branch in the (y, z)-plane becomes more 
complex because the fundamental V\ develops a critical 




FIG. 3: Contours of \v ± \ for (a, 7) = (1,2) at R = 50171 
(solid) and 12637 (dashed) stretched by R 1 ^ 3 factors along 
curves normal to ito-contours to match \vi\ contours at R = 
3079 (dash-dot). The (almost overlapping) black and yellow 
solid curves show uo(y, z) = at the 3 R's. 



layer, as discussed hereafter. 

Figure [5] illustrates the structure of the lower branch 
steady state. The streaky flow u n (y,z) and the rolls 
(0, Vq, wq) remain large scale and their structure becomes 
independent of R; vq has a simple updraft at z = and 
downdraft at z = ±ir/j that sustain the z modulation 
of Uo(y, z) (recall that uq = ±1 at y = ±1 in PCF). But 
the fundamental mode v\ concentrates about the critical 
layer uo(y, z) = c (c = for these states in PCF). Critical 
layers are well-known in the context of the 2D, linear the- 
ory of shear flows [l]| . Here the critical layer is a surface 
in 3D space and it is nonlincarly coupled to the 0-mode 
Vo(y,z). When the higher harmonics become negligible 
and I ^ol, I Wo I <C |unL the equation for the fundamental 
o [11, 



mode simplifies to 



14 



[ia(u - c)vi + (vi ■ Vuo)e x ] e ie = 

- V( Pl e w ) + R-'V 2 { Vl e ie ) , (3) 

with V • (vie l6 ^j = 0. For high R, the solutions develop 
an i? -1 / 3 critical layer in the neighborhood of Uo(y, z) — 
c = that results from the balance between a(uQ — c) = 
a(r - r c ) ■ Vit + 0(\r — r c \ 2 ) and i? _1 V 2 , so if 5 is the 
critical layer thickness, we must have aS\ Vuo| ~ R~ 1 5~ 2 
and 8 ~ (cv| Vuol-R)" 1 ^ 3 near Uo(y, z) — c = 0. Figure [3] 
confirms that critical layer scaling for the lower branch 
steady state in PCF. 

The nonlinear coupling 0, [ll[ between the fundamen- 
tal Vi, with its critical layer structure, and the rolls 
(vo,Wq) provides a challenge for the development of a 
full asymptotic theory of the lower branch states that 
would be able to predict the amplitude scaling of the 
fundamental mode. If v\ remained a large scale struc- 
ture, its amplitude would have to scale like R~ l in order 
for its nonlinear self-interaction to balance the viscous 
diffusion of the strcamwise rolls vq jTl , 14 1 . The de- 
velopment of a critical layer scale complicates the anal- 
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ysis and different norms and components have different 
scalings. Nonetheless, an asymptotic theory appears fea- 
sible and the present numerical data is clear and its im- 
plications are significant: the lower branch states tend 
to a relatively simple but non-trivial quasi-2D singular 
asymptotic state as R — > oo that is not a solution of the 
Euler equation (eqn. ([1]) with = 0), and that is not 
the laminar flow v = ye x either. So the lower branch 
states do not bifurcate from the laminar flow, not even 
at R = oo. The data presented is for (a, 7) = (1,2) 
however identical features hold for other (a, 7) values. 

Turning now to a stability analysis of the lower branch 
coherent states we find that these states are distinguished 
not only by their asymptotic structure but also by their 
stability characteristics. Our eigenmode analysis of the 
3D lower branch steady state in plane Couette flow, up 
to R = 12000, show that they have a single, real unstable 
eigenvalue shown in figurc[4]for (a, 7) = (1.14,2.5). This 
state is most unstable at R w 342 then the unstable 
eigenvalue steadily decreases approximately as R~ 0A8 for 
larger R. Furthermore, the corresponding eigenfunction 
is in the same shift-reflect and shift-rotate symmetries 
@, eqns. (24), (26)] as the lower branch state. This is 
not true for the upper branch states which develop new 
bifurcations and unstable modes as R increases. 

These stability results were obtained using both a di- 
rect calculation of the eigenvalues of the full Jacobian in 
the doubly symmetric subspace of the lower branch state 
with an ellipsoidal truncation of the Fourier-Chebyshcv 
representation Q , and an iterative calculation in the full 
space using the Arnoldi algorithm and the ChannelFlow 
code with cubic truncation [3]. The leading unstable 
and least stable eigenvalues matched to 5 or 6 signifi- 
cant digits. We have also investigated subharmonic insta- 
bilities through numerical simulations in a double-sized 
box with fundamental wavenumbers a/2 and 7/2. Sam- 
ple simulations with 'random' perturbations did not re- 
veal further instabilities, however a more systematic ap- 
proach using the Arnoldi algorithm revealed a weakly 
unstable subharmonic in x. For R = 1000 and (a, 7) = 
(1.14,2.5), the fundamental instability shown in fig. 0] 
has growth rate 0.03681 while the subharmonic instabil- 
ity has growth rate 0.005248 ± i 0.02245. The analysis of 
this subharmonic mode is left for future study. 

Thus, in the one-period domain with fundamental 
wavenumbers (a, 7), the lower branch state is an unsta- 
ble equilibrium with a ID unstable manifold. Therefore 
its stable manifold splits the phase space into two parts, 
at least locally. The evolution of disturbances in the one- 
period domain, starting on the ID unstable manifold of 
the lower branch on either side of the stable manifold 
is illustrated in fig. [5] These numerical simulations were 
performed using ChannelFlow in the full phase space and 
show the time evolutions in the energy input-energy dis- 
sipation plane, both normalized by their laminar values. 
For plane Couette flow, the normalized energy input rate 
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FIG. 4: The single unstable eigenvalue of the lower branch 
state (a, 7) = (1.14, 2.5) as function of R. Asympotic scaling 
is w O(R~ 0A8 ). There is an extra complex conjugate pair 
near the onset R an ~ 218. 



is equal to the normalized drag, that is, the drag at the 
walls normalized by their laminar value. Perturbations 
starting on one side of the stable manifold gently de- 
cay back to the linearly stable laminar flow v = ye x 
while perturbations on the other side of the stable man- 
ifold shoot to a turbulent state. Figure O also shows the 
upper branch sister of the lower branch state which, as 
stated earlier, is located in phase space much closer to the 
'turbulent' state. The decay of perturbed lower branch 
states back to the laminar flow follows a standard two- 
step evolution. First, the fundamental mode Vi, with 
its critical layer structure, disappears and the flow re- 
laxes to an ^-independent state that consists of stream- 
wise rolls (0,vq,wo) and streaks u (y,z) and slowly de- 
cays back to the laminar flow on a long viscous time 
scale. Perturbations that shoot to a turbulent state fol- 
low a much more rapid 'breakdown' with high dissipation 
rate (about 13 on fig. then settle to a turbulent state 
with energy input and dissipation rates of about 4.4 (for 
(a,j,R) = (1,2,1000)). 

These results suggest that the lower branch stable 
manifold is the boundary separating the basin of attrac- 
tion of the laminar state from that of the turbulent state 
and therefore that they may be the key states control- 
ling transition to turbulence. Our results have focused 
on symmetric steady states in plane Couette flow but 
there is evidence of a similar role for lower branch trav- 
eling waves in plane Poiseuille flow [B| and pipe flow 17 1 . 
Recent work by Viswanath [l8| complements our work 
by showing that perturbations of the laminar flow in the 
form of streamwise rolls of the right threshold amplitude 
+ small 3D noise do get attracted to a lower branch state 
before shooting to turbulence. We expect the symmet- 
ric lower branch state to play a key role for transition 
in plane Couette flow but there exist other asymmet- 
ric lower branch traveling wave states as well as peri- 



odic orbits [18|, |19|, |2fJ , each of which may play a simi- 



lar 'transition-backbone' role, locally in phase space. We 
conjecture that the permanent states (steady states, trav- 
eling waves and periodic orbits) most relevant to transi- 
tion to turbulence will contain streamwise rolls. It 
may be possible to trigger transition with smaller dis- 
turbances but we suspect that such disturbances would 
necessarily lead to the formation of i?" 1 rolls and ap- 
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FIG. 5: Energy input /dissipation rate starting near the lower 
branch fixed point (a, 7, R) = (1, 2, 1000) on its unstable man- 
ifold. In one direction, the flow goes to turbulence while in the 
other direction it relaminarizes. The dot spacing is At = 5. 
The blue marker at (1,1) is the laminar flow, green marker 
at (1.35,1.35) is the lower branch state and the red marker at 
(3.89,3.89) is its upper branch sister. 



proach toward a lower branch state with _R _1 rolls prior 
to transition along the unstable manifold of the lower 
branch state. 

The extreme low dimensionality of the lower branch 
unstable manifold suggests a new approach to turbulence 
control. Turbulence control strategies roughly fall into 
2 categories: either prevent nonlinear breakdown of the 
linearly stable laminar flow, or push the fully nonlinear 
turbulent flow back to laminar. A new strategy might be 
to put the flow on the lower branch equilibria and keep it 
there by controlling its very few unstable modes. There 
is a small drag penalty to do so since lower branch states 
have a net drag that is 30 to 40% higher than the lami- 
nar state as R — ► 00 but that is a massive drag reduction 
compared to the turbulent state. This control strategy is 
related to, but quite distinct from the strategies proposed 
Streaks are used in to efficiently de- 
revent the linear 



2l| and [25 



form the laminar base flow in order to 
instability of boundary layer flow. In [2l| strategies are 
considered to push the turbulent flow onto the laminar 
side of the stable manifold of a lower branch unstable 
periodic solution in order to relaminarizc the flow. The 
current proposal is to put the flow on the unstable lower 
branch equilibrium and keep it there by controlling its 
single unstable cigenmode. 

JW and FW were partially supported by NSF grant 
DMS-0204636. We thank Divakar Viswanath and Pre- 
drag Cvitanovic for helpful discussions. 



t Electronic address: gibson@cns.physics.gatech.edu 
* Electronic address: waleffe@math.wisc.edu 
[1] B. Hof, A. Juel, and T. Mullin. Scaling of the turbu- 
lence transition threshold in a pipe. Phys. Rev. Lett., 
91:244502, 2003. 
[2] R. Fitzgerald. New experiments set the scale for the onset 
of turbulence in pipe flow. Physics Today, 57(2):21-23, 
2004. 

[3] F. Waleffe. Three-dimensional coherent states in plane 
shear flows. Phys. Rev. Lett., 81:4140-4148, 1998. 

[4] F. Waleffe. Exact coherent structures in channel flow. J. 
Fluid Mech., 435:93-102, 2001. 

[5] T. Itano and S. Toh. The dynamics of bursting process in 
wall turbulence. J. Phys. Soc. Japan, 70:703-716, 2001. 

[6] F. Waleffe. Homotopy of exact coherent structures in 
plane shear flows. Phys. Fluids, 15:1517-1543, 2003. 

[7] H. Faisst and B. Eckhardt. Traveling waves in pipe flow. 
Phys. Rev. Lett, 91:224502, 2003. 

[8] H. Wedin and R.R. Kerswell. Exact coherent structures 
in pipe flow. J. Fluid Mech., 508:333-371, 2004. 

[9] J. Jimenez, G. Kawahara, M.P. Simens, and M. Nagata. 
Characterization of near- wall turbulence in terms of equi- 
librium and 'bursting' solutions. Phys. Fluids, 17:015105 
(16pp.), 2005. 

[10] B. Hof, C.W.H. van Doorne, J. Westerweel, F.T.M. 
Nieuwstadt, H. Faisst, B. Eckhardt, H. Wedin, R.R. Ker- 
swell, and F. Waleffe. Experimental observation of non- 
linear traveling waves in turbulent pipe flow. Science, 
305(5690):1594-1598, 2004. 

[11] F. Waleffe. On a self-sustaining process in shear flows. 
Phys. Fluids, 9:883-900, 1997. 

[12] G. Kreiss, A. Lundbladh, and D. S. Henningson. Bounds 
for threshold amplitudes in subcritical shear flows. Jour- 
nal of Fluid Mechanics, 270:175-198, July 1994. 

[13] S.A. Maslowe. Critical layers in shear flows. Ann. Rev. 
Fluid Mech., 18:405-432, 1986. 

[14] F. Waleffe. Hydrodynamic stability and turbulence: Be- 
yond transients to a self-sustaining process. Stud. Applied 
Math., 95:319-343, 1995. 

[15] D.J. Benney. The evolution of disturbances in shear flows 
at high Reynolds numbers. Stud. Appl. Math., 70:1-19, 
1984. 

[16] J. F. Gibson. Channelflow: a spectral Navier-Stokes sim- 
ulator in C++. Technical report, Georgia Institute of 
Technology, 2006. http://www.channelflow.org. 

[17] R.R. Kerswell and O.R. Tutty. Recurrence of traveling 
waves in transitional pipe flow. J. Fluid Mechanics (sub- 
mitted), arxiv.org/physics/0611009, 2007. 

[18] D. Viswanath. The dynamics of transition to turbulence 
in plane Couette flow, arxiv.org/physics/0701337, 2007. 

[19] D. Viswanath. Recurrent motions within plane Couette 
turbulence. J. Fluid Mech., 2007 (to appear). 

[20] G. Kawahara and S. Kida. Periodic motion embedded in 
Plane Couette turbulence: regeneration cycle and burst. 
J. Fluid Mech., 449:291-300, 2001. 

[21] G. Kawahara. Laminarization of minimal plane Couette 
flow: Going beyond the basin of attraction of turbulence. 
Physics of Fluids, 17(4):041702, 2005. 

[22] J.H.M. Fransson, A. Talamelli, L. Brandt, and C. Cossu. 
Delaying transition to turbulence by a passive mecha- 
nism. Physical Review Letters, 96(6):064501, 2006. 



* Electronic address: |wang@math.wis c.edu| 



